function diff = model_shock_SA( theta, mean, param, aux)

    param.lambda = exp(theta);
	param = model_SA(param, aux);    
    diff = (param.mean_x - mean);
    
end